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Abstract 

The static stress needed to depin a 2D edge dislocation, the lower dynamic stress needed to keep 
it moving, its velocity and displacement vector profile are calculated from first principles. We use 
a simplified discrete model whose far field distortion tensor decays algebraically with distance as in 
the usual elasticity. Dislocation depinning in the strongly overdamped case (including the effect of 
fluctuations) is analytically described. N parallel edge dislocations whose average inter-dislocation 
distance divided by the Burgers vector of a single dislocation is L ^ 1 can depin a given one if 
N = 0{L). Then a limiting dislocation density can be defined and calculated in simple cases. 

PACS numbers: 61.72.Bb, 5.45.-a, 45.05.+X, 82.40.Bj 



1 



In many fields, genuinely microscopic phenomena affect macroscopic behavior in a way 
that is difficult to quantify precisely. Topical cases are the motions of dislocations |l|, |2| , 
cracks [sl, vortices in Josephson arrays 4] or other defects subject to pinning due to the un- 
derlying crystal microstructure. Emerging behavior due to motion and interaction of defects 
might explain common but poorly understood phenomena such as friction jsj] . Macroscopic 
theories consider the continuum mechanics of these solids subject to forces due to the defects 
and additional equations for the densities of defects and properties of their motion Q]. The 
latter are usually postulated by phenomenological considerations. An important problem is 
to derive a consistent macroscopic description taking into account the microstructure. 

Here we tackle a simplified problem containing all the ingredients of the previous de- 
scription: the pinningand motion of edge dislocations. Firstly, we study a two-dimensional 
(2D) discrete model ['i^ describing the damped displacement of atoms subject to the field 
generated by a 2D edge dislocation and a constant applied shear stress of strength F. If 
|F| < Fes {Fes is related to the static Peierls stress), the stable displacement field is station- 
ary, whereas the dislocation core and its surrounding displacement field move if \F\ > Fes- 
A crucial observation is that there exists a stable uniformly moving dislocation with both 
core and far field advancing at the same constant velocity. This suggests that a moving edge 
dislocation is a traveling wave of the discrete model. 

Our self-consistent calculation ^ based on this picture predicts the following magnitudes: 
(i) the critical static stress needed to depin a static dislocation, (ii) the dynamic stress Fed < 
Fes below which a moving dislocation stops (in the strongly overdamped case. Fed = Fes), 
and (iii) the dislocation velocity as a function of applied stress. The latter information has 
to be taken from experiments in the standard treatment |j9|. Then a macroscopic quantity, 
the dislocation velocity, is obtained from analysis of a microscopic model. 

Secondly, we consider a distribution of many parallel edge dislocations separated by 
macroscopic distances comprising many lattice periods. A dislocation cannot move under the 
influence of other dislocations far away unless the latter have finite density (there are N such 
dislocations and the average distance between them is L with = 0{L) as L ^ oo). Under 
the influence of such a distribution, one dislocation may be pinned or move depending on the 
dislocation density. The latter is calculated at the critical stress in a simple configuration. 

A simplified discrete model of edge dislocations. Consider an infinite three dimensional 
cubic lattice with symmetry axes x, y, z. We insert an extra half plane of atoms parallel to 
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the plane yz. The border of this extra half plane is a line (parallel to the z axis), which is 
called an edge dislocation in the crystal. The Burgers vector of the dislocation points in the 
X direction and the plane xz is the glide plane of the dislocation (see 0]). If we apply an 
external stress a, the dislocation moves on its glide plane and in the direction of its Burgers 
vector in response to just one component of a: the stress r resolved on the glide plane in 
the glide direction 3, 0|. All the sections of the lattice by planes parallel to xy look alike. 
Thus, we reduce the problem to a 2D lattice with an extra half line of atoms, see Fig. 22 
of Ref. y. Assuming that glide is only possible in the x-direction, the dynamics of an edge 
dislocation in a 2D lattice can be described by Q]: 

d Ui j dui j 

^ ~w ~dr ^ ~ 

+ A[sm{uij+i-Uij) + sm{uij-i-Uij)]. (1) 

The lattice is a collection of chains in the x direction with elastic interaction between nearest 
neighbors within the same chain and sinusoidal interaction between chains. Uij/ {2tt) is the 
dimensionless displacement of atom (^, j) in the x direction measured in units of the Burgers 
vector length b. A > measures the relative strengths of the nonlinear forces exerted 
by atoms on different planes y = k (constant) and the linear forces exerted within any 
plane y = k. The dimensionless parameter A also determines the width of the dislocation 
core (~ I/a/A). Lastly, the time unit is the ratio between the friction coefficient and the 
spring constant in the x direction. Then m is the dimensionless ratio between the atomic 
mass times the spring coefficient and the square of the friction coefficient. In dislocation 
dynamics, an important case is that of overdamped dynamics, m = Eq. (Q) can 



be generalized to a vector model having a displacement vector 



'uij,Vij) and a continuum 
Such model has among 



limit yielding the 2D Navier equations with cubic symmetry 
its solutions edge dislocations with Burgers vectors in the x or y directions gliding in the 
direction thereof (which does not have to be assumed as in the present simple model). This 
model can also be solved using our methods at the expense of technical complications and 
high computational cost. 

In this geometry, the far field of a static 2D edge dislocation is approximately given by 
the corresponding continuum elastic displacement Uij = u{x,y) with x = ei, y = ej (where 
e = b/ L <^ 1, i, j are large and L is the appropriate mesoscopic length). Then the stationary 
solutions of Eq. ((H) satisfy the equations of anisotropic linear elasticity, u^x + ^Uyy = u^x + 



Uyy = 0, (Y = y/ ^/A) far away from singularities and jumps The solution corresponding 
to the edge dislocation is the polar angle 9{x,Y) G [0,27r), measured from the positive x 
axis. Continuum approximations break down near the dislocation core, which should be 
described by the discrete model Q]. The advantadge of Eq. ((T)) compared to other 2D 
generalizations [13] of the Frenkel-Kontorova model is that it yields the correct decay for 
strains and stresses: r""*^ as = + ^ oo, instead of exponential decay. 

Overdamped dynamics and static Peierls stress. We shall now study the structure of 
a static edge dislocation of Eq. ((H), the critical stress needed to set it in motion and its 
subsequent speed. We solve numerically Eq. (jT)) with m = on a large lattice < N 

using Neumann boundary conditions (NBC) corresponding to applying a shear stress of 
strength F in the x direction. The (far field) continuum elastic displacement for a static 
2D edge dislocation subject to such a shear stress is 6(x,y /^/A) + Fy. Then the NBC are 

U±{N+l),j - U±N,j = ±[^±{N+l),j - ^±N,j] and Ui^±(N+l) - Ui,±N = ±[^t±iN+l) - ^t±N] ± ^ , 

where = 9{i,j/\fA) with 6^(0,0) = tx/2. If F = and the initial condition is the elastic 
far field Of^, the system relaxes to a stationary configuration Uj The dislocation is expected 
to remain stationary for F 7^ unless \F\ is larger than a critical value Fcs{A), related to 

3| • Nonlinear stability of the stationary edge dislocation 
14 To test this picture, we solve numerically Eq. ([Q) in a 
large lattice, using NBC and the static dislocation obtained for F = as initial data. For 
large times and |F| small, the system relaxes to a steady configuration Ui^j which provides 
the structure of the core, see Fig. H When |F| is large enough, the dislocation is observed 
to glide in the x direction: to the right if F > 0, and to the left if F < 0. 



the so called static Peierls stress 
for |F| < Fes was proven in Ref. 



To calculate F^^, we extend the depinning calculations of Ref. Il5| to 2D systems. We 
redefine = Uij + Fj, insert f/jj = Uij{F, A) + Vij{t) in Eq. ^ with m = and expand 
the resulting equation in powers of Vij, about the stationary state Uij{A,F) up to cubic 
terms. Subscripts in the resulting equation can be numbered with a single one starting from 
the point i = j = -N: Uij = Uk and Vij = Vk, k = i + + N) {2N + 1) for z, j = -N, ...,N. 
The resulting equation can be written formally as dw/dt = A4{F)v + i3(v, v; F), where the 
vector V has components Vk- The linear stability of the stationary state Uk{A,F) depends 
on the eigenvalues of the matrix Ai{F). These eigenvalues are all real negative for |F| < Fes 
whereas one of them vanishes at |F| = F^s- This criterion allows us to numerically determine 
Fes as a function of A; see Fig. E^a). Notice that the critical stress increases with A. Thus 
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narrow core dislocations {A large) are harder to move. 

Dislocation velocity. Let us assume that F > Fes (the case F < —Fes is similar). Then 



{F-Fe)j + 
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tjrvj- (plus terms that decay exponentially fast in time). The procedure 
lyfor discrete ID systems yields the amplitude equation d(f)/dt = a + 



sketched in Refs. 

Here a = l-Ar(F,,) (F-F,,)/(l-r), = A [cos(f/,,,-+i-f/,,,-+F,)-cos(f/,,,-f/i,,-_i+F,)] 
and f3 = 1 ■ i3(r, r; Fes) / (1 ■ r) , 1 and r are the left and right eigenvectors of the matrix Ai (Fes) 
corresponding to its zero eigenvalue (its largest one). From the amplitude equation, the 
approximate dislocation velocity is ^15^: c ~ a/o^/tt = 0{\F—Fcs\^). Numerically measured 
and theoretically predicted dislocation velocities are compared in Fig. El^b). Calculations in 
lattices of different sizes yield similar results. 

How do we calculate numerically the dislocation velocity? This is an important point 
for using the calculated dislocation velocity as a function of stress in mesoscopic theories 
and a few comments are in order. If we solve numerically Eq. with static NBC for 
\F\ > Fes, the velocity of the dislocation increases as it moves towards the boundary. The 
dislocation accelerates because we are using the far field of a steady dislocation as boundary 
condition, instead of the (more sensible) far field of a moving dislocation. However, the 
latter is in principle unknown because we do not know the dislocation speed. We will assume 
nevertheless that the dislocation moves at constant speed c once it starts moving, as it would 
in a stressed infinite system. Then the correct dislocation far field is 6(i — ct,j/^/A) + Fj. 
With this far field in the NBC, Eq. has traveling wave solutions Uij{t) whose velocity 
can be calculated self-consistently. How? By an iterative procedure that adopts as initial 
trial velocity that of a dislocation subject to static NBC as it starts moving. Near threshold, 
step-like profiles are observed (see Fig.O)), that become smoother as F increases. The profiles 
have been calculated by following the trajectories of points with the same value of i and 
different values of j for F > Fcs{A), according to the formula Uij(t) = u{(,j), ( = i — ct. 
Notice that the wave front profiles are kinks for j < and antikinks for j > 0. 

Influence of fluctuations. The original discrete model contains both damping and fluctu- 



ation tenns B. V.uctuat.o. te™. a. app.ecab.e o„,v .ea. and cont.bute a. add.tive 

white noise term to the amplitude equation. Due to this term, there is a small probability 
for the dislocation to move even if \F\ < F^s and m = 0. The resulting average velocity can 
be estimated by observing that the potential of the corresponding Fokker-Planck equation 

3 

is cubic and it has a small barrier of height proportional to {F^s — |-F|)2 . Then the exponen- 
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tially small velocity of the dislocation under the critical stress is the reciprocal of the mean 
escape time from the barrier Provided {Fes — ^ D (where D measures the noise 
strength), we have — Inc oc {F^s — \F\)2/D. 

Inertia and dynamic Peierls stress. Inertia changes the previous picture of dislocation 
motion in one important aspect: the dislocations keep moving for an interval of stresses 
below the static Peierls stress, F^ < |F| < Fes- On this stress interval, stable solutions 
representing static and moving dislocations coexist: to depin a static dislocation, we need 
\F\ > Fes- However if \F\ decreases below F^s, a moving dislocation keeps moving until 
|F| < Fed, see Fig. EI Thus Fed 

Our theory therefore yields the static and the dynamic Peierls stresses and the velocity of a 
dislocation. 

Interaction between edge dislocations. Let us assume that there are static edge dislo- 
cations at the points {xn,yn) parallel to one dislocation at {xo,yo), and that all dislocations 
are separated from each other by distances of order L ^ 1 (measured in units of the Burg- 
ers vector length, b). We want to analyze whether the collective influence of the N distant 
dislocations can move that at {xo,yo). This problem is similar to that of deriving a reduced 
dynamics for the centers of 2D vortices of Ginzburg-Landau equations subject to their mu- 
tual influence 3|- In the case of dislocations, the existence of a pinning threshold implies 
that the reduced dynamics is that of a single dislocation subject to the mean field created 
by the others. We thus have a reduced field dynamics, not particle dynamics as in the case 
of the Ginzburg-Landau vortices. 

Displacement vectors pose the problem of defining branch cuts in the continuum (elastic) 
limit, which leads us to consider instead the distortion tensor as a primary quantity [6]. For 
our discrete system, the distortion tensor has nonzero components w^^j = Uj+ij — Uij and 
Wj-j = sin(Mjj_|_i — Uij) that become e (du/dx) and e (du/dy), respectively, in the continuum 
limit e = 1/L 0, x—xq = ei, y—yo = ej finite. In the continuum limit, the distortion tensor 
of an edge dislocation centered at the origin has nonzero components w'-^-* = —e\^y/(Ax'^ + 
y"^) and w*^^-* = e\fAx/{Ax'^ + y'^). If we have edge dislocations at {xn,yn), 1 < n < N, 
far from one at {xo,yQ), the distortion tensor is sum of individual contributions. Then the 
far field distortion tensor seen by the dislocation at (xo,?/o) is: 



(1) , TP (2) ^/^^ , IP fn\ 
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F, = ey ~ + .... (4) 

The dislocation at (xo,?/o) moves if F2 > Fcs{A). This cannot be achieved as e — unless 
= 0(l/e). Then the sums in Eqs. (jH)) and become integrals. We define a static 
dislocation density p{x, y) as the limit of A^~^ Hn=i "^(^ ~ ^n)S{y — yn) as N ^ 00. Then 

A{xo - xy + (yo - y) 



00 roo 



00 •! —00 



' ' A{x^-xf + {y^-yf ^ ' 



—00 ^ —00 



where a = eN is the ratio of the total Burgers vector to the mesoscopic length measuring 
average interdislocation distance. As an example, let us assume that yo = yn = 0. Then 
p = p{x)6{y) and Fi = 0. Let us assume that the dislocations are constrained by two 
obstacles at xo = ±/ and subject to the same critical stress. Then the critical dislocation 
density is p{x) = [1 — {^/AFcs/(y) x]/['n-y/P — x^], provided a > y/AFcsl (cf. sl page 127). 

In conclusion, edge dislocations can be characterized as traveling waves of discrete mod- 
els. The dislocation far field moves at a constant velocity equal to that of the dislocation 
core. Static and dynamic Peierls stresses and the dislocation velocity as a function of ap- 
plied stress can be found numerically (or analytically near critical stress in the overdamped 
case) and adopted as the basis of a mesoscopic theory [l^. We have also shown that the 
interaction between distant edge dislocations can be described in terms of a continuous dis- 
location density. This field-theoretical reduced description greatly contrasts with the case 
of interacting point vortices, which is completely described by the particle dynamics of the 
vortex centers Extension to fully vectorial models and to other types of dislocation 

should follow along similar lines. 
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FIG. 1: Displacement field profile for the stationary edge dislocation with A = 1 and N = 50. 




FIG. 2: (a) Static (squares, m = 0) and dynamic (asterisks, m = 0.5) critical stresses Fes and Fed 
versus A. (b) Theoretical (solid line, m = 0) and numerical (squares, m = 0; asterisks, m = 0.5) 
dislocation velocity vs. F {A = 1, N = 25). 
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